Global fits of the Minimal Universal Extra Dimensions scenario 
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In theories with Universal Extra-Dimensions (UED), the 71 particle, first excited state of the hy- 
percharge gauge boson, provides an excellent Dark Matter (DM) candidate. Here we use a modified 
version of the SuperBayeS code to perform a Bayesian analysis of the minimal UED scenario, in 
order to assess its detectability at accelerators and with DM experiments. We derive in particular 
the most probable range of mass and scattering cross sections off nucleons, keeping into account 
cosmological and electroweak precision constraints. The consequences for the detectability of the 
71 with direct and indirect experiments are dramatic. The spin-independent cross section proba- 
bility distribution peaks at ~ I0 -11 pb, i.e. below the sensitivity of ton-scale experiments. The 
spin-dependent cross-section drives the predicted neutrino flux from the center of the Sun below the 
reach of present and upcoming experiments. The only strategy that remains open appears to be 
direct detection with ton-scale experiments sensitive to spin-dependent cross-sections. On the other 
hand, the LHC with I fb^ 1 of data should be able to probe the current best-fit UED parameters. 
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I. INTRODUCTION. 

Dark matter (DM) studies are often carried out in 
the framework of Supersymmetric (SUSY) theories, but 
there are many alternative extensions of the Standard 
model of particle physics that lead to viable DM can- 
didates. Among them, models with universal extra di- 
mensions (UED), in which all Standard Model (SM) par- 
ticles can propagate in the bulk of one or more com- 
pactified flat extra dimensions [l[ , have received consider- 
able attention, and they have been studied in relation to 
collider phenomenology, indirect low-energy constraints, 
cosmology, and dark matter (see Refs. 0, 3 for recent 
reviews). In the simplest and most popular version, there 
is a single extra dimension compactified on an interval, 
S\/Z2- Each SM particle has a whole tower of Kaluza- 
Klein (KK) modes, labelled by an integer n, called KK 
number, which is nothing but the number of quantum 
units of momentum which the SM particle carries along 
the extra dimension. 

One of the peculiar features of UED theories is the 
conservation of the KK number at tree level, which is a 
simple consequence of momentum conservation along the 
extra dimension. This implies that the lightest KK parity 
odd particle (LKP) is stable over cosmological timescales, 
and being cold and neutral, it provides a suitable WIMP 
candidate 0-[3. 

In this paper we shall concentrate on the Minimal Uni- 
versal Extra Dimensions (MUED) discussed in Ref. @, 
where there are only two parameters in addition to the 
Higgs mass (rrih), namely the size of extra dimension (R) 
and the cutoff scale of the theory, A (see Refs. (l0l - fr7| for 
the case of two universal extra dimensions) . Often AR is 
used instead of the cutoff scale itself, and we follow the 
same convention. 



The LKP in MUED turns out to be the level- 1 KK 
partner (71) of the SM photon, which must therefore be 
stable by virtue of the conservation of KK parity. In 
the limit where > v, where v is the vacuum ex- 

pectation value of the Higgs, the Weinberg angle for KK 
gauge boson is negligible, and the KK gauge bosons are in 
fact almost weak eigenstate. Therefore the KK photon 
is almost the hypercharge KK gauge boson (7„ ps B n ) 
and the KK Z is almost the neutral weak eigenstate of 
SU(2) W (Z n » Wl). 

While we do not know much about the DM particle, 
there are many direct, indirect and accelerator searches 
current undergoing, with the aim of identifying them 
PH-Hol. KK DM provides a valid alternative to the 
widely discussed Supersymmetric DM [2l[ , and it is often 
adopted as a case study scenario when trying to assess 
the capability of experimental strategies to discriminate 
among various DM candidates (e.g. [22|). 

In this paper we perform a Bayesian analysis of the 
MUED scenario, in order to assess its detectability at 
accelerators and with DM experiments. We derive in par- 
ticular the most probable range of mass and scattering 
cross sections off nucleons, keeping into account cosmo- 
logical and electroweak precision constraints. As we shall 
see, this has dramatic implications for the detectability 
of KK DM. 

The paper is organized as follows: in section [TT| we dis- 
cuss the theoretical framework of the MUED scenario. 
In section IIIII we provide some details on our statistical 
tools, including a discussion of the priors adopted in the 
Bayesian analysis of the MUED parameter space. In Sec. 
IIVI we present the results and in|V]we discuss their con- 
sequences and conclude. 
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II. THEORETICAL FRAMEWORK 
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In the MUED, the vanishing boundary conditions are 
assumed for all KK particles at the cutoff scale (i.e., all 
KK particles at level-n are degenerate (m n = n/R) at the 
cutoff scale), and therefore the mass spectrum at elec- 
troweak (EW) scale is completely determined by 

RG evolution between and A [4[ (there is also a con- 
tribution from EW symmetry breaking, which is small 
except for top quark). Since the estimated cutoff scale 
is not too far away from the resulting mass spec- 

trum is somewhat degenerate due to short RG running. 
As expected from RG running, the masses of the KK 
particles depend on how strongly they interact, therefore 
strongly interacting KK particles get larger corrections 
than weakly interacting particles. In MUED, the KK 
gluon is the heaviest particle, followed by KK quarks, 
KK Z/W and KK leptons. 

Due to KK-parity, contributions to clcctroweak observ- 
ables do not appear at tree-level and this allows KK par- 
ticles to be light enough so that they can be produced 
at current collider experiments. This has been studied 
in Refs. 0, [23| and revisited more recently in Ref. [24| 
including subleading contributions as well as two loop 
corrections to the SM p parameter. A lower bound on 
R^ 1 from those oblique corrections is ^600 GeV at 90% 
C.L. with a Higgs mass of 115 GeV, which is the LEP 
limit. However this constraint is significantly relaxed 
with increasing Higgs mass, allowing for a compactifica- 
tion scale as low as 300 GeV (other indirect low-energy 
constraints are comparable or weaker. See 0, Q and ref- 
erences therein.). 

In the rest of this section, we briefly review the calcu- 
lation of the relic density. Since KK particles in MUED 
are somewhat degenerate, it is important to include coan- 
nihilation effects. The generalization of the relic den- 
sity calculation including coannihilations is straightfor- 
ward [j| H^] . Assume that the particles \i are labelled 
according to their masses, so that m, < nij when i < j. 
The number densities rii of the various species Xi obey a 
set of Boltzmann equations. It can be shown that under 
reasonable assumptions (25| . the ultimate relic density 
n x of the lightest species xi (after all heavier particles 
Xi have decayed into it) obeys the following simple Boltz- 
mann equation 



dn x 
~dt 



(<r eS v)(nl -nL) 
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where H is the Hubble parameter, v is the relative veloc- 
ity between the two incoming particles, n eq is the equi- 
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Here <7y = cr(xiXj ~^ SM) are the various pair annihila- 
tion cross sections into final states with SM particles, gi 
is the number of internal degrees of freedom of particle 
Xi and n x = ^ t n, is the density of xi we want to 
calculate. 

By solving the Boltzmann equation analytically with 
appropriate approximations [H, HH , the abundance of the 
lightest species xi ls given by 
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where the Planck mass scale is Mpi = 1.22 x 10 19 GeV 
and is the total number of effectively masslcss degrees 
of freedom at temperature T: 

9*{T)= 9i + l 9i ■ 

i— bosons i— fermions 
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where the constant c is determined empirically by com- 
paring to numerical solutions of the Boltzmann equation 
and here we take c = ^ as usual. a c s and 6 e ff are the first 
two terms in the velocity expansion of a e e 



a e ff(x) v = a cB (x) + b cS (x) v 2 + 0(v A ) . 
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where Oy and bij are obtained from aijV = a^- + bijV + 
0(v A ) and v is the relative velocity between the two an- 
nihilating particles in the initial state. Considering rela- 
tivistic corrections to the above treatment results in an 
additional sublcading term which can be accounted for 
by the simple replacement 

b^b-^a, (13) 

in the above formulas. For our calculation of the relic 
density, wc use the cross sections given in Refs. [1, HH, H3 ■ 
For resonance effect, which we do not include, see Refs. 
[H-[30j]. In the MUED, coannihilation with SU(2) W - 
singlct KK leptons (ei) is important since it is the next- 
to-lightest KK particle, and mei ~" 71 - 0.01. 



III. STATISTICAL ANALYSIS 

The free parameters of the model are the SM Higgs 
mass, rrih, the inverse radius of the UED, i?" 1 , and the 
cutoff scale A. For numerical reasons, we work with the 
following MUED parameters: 

*= {m h ,R-\AR}. (14) 

In particular, we adopt as a free parameter the number 
of KK levels AR rather than A itself. In our scan, AR is 
considered as a real- valued variable, but we the round it 
to the nearest integer value when computing the observ- 
able quantities. We also include in our scan as nuisance 
parameters the relevant SM parameter set 

$={M t , m b {m b Y^, aem(M z p, a s (Af z )^} , 

(15) 

where M t is the pole top quark mass, while the other 
three parameters (the bottom mass, the electromagnetic 
and the strong coupling constants) are all evaluated in 
the MS scheme at the indicated scales. 

We denote by 8 = {^f, <&} the vector of parameters en- 
tering the analysis, and by d the available data (described 
below). Bayes theorem reads 

Ptm = £6™, ,16, 

where P(Q\d) is the posterior distribution on the param- 
eters (after the observations have been taken into ac- 
count), P(d\Q) = C(<9) is the likelihood function (when 
considered as a function of 9 for fixed data d) and P{&) 
is the prior distribution, which encompasses our state of 
knowledge about the value of the parameters before we 
have seen the data. Finally, the quantity in the denomi- 
nator of Eq. (|T5j) is the Bayesian evidence (or model like- 
lihood), a normalizing constant which does not depend 
on and which can be neglected when one is interested 
in parameter inference. Together with the model we must 
specify the priors for the parameters, which enter Baycs' 



theorem, Eq. ([i~6|) . As in any good Bayesian analysis it 
is important to asses the relevance of prior choices, we 
perform our scan using two different priors: 

• Flat prior: a uniform prior over the ranges 
10 GeV < m h < 3 TcV, 280 GeV < R- 1 < 3 
TeV and 1 < AR < 100. 

• Log prior: a uniform prior over the ranges 1 < 
log(m/,/GeV) < 3.5, 2.4 < log(i?-7GcV) < 3.5 
and 1 < AR < 100. 

The lower bound on comes from considering current 
collider limits from trilcpton search at the Tevatron, giv- 
ing R- 1 > 280 GeV at 95% C.L. with 100 pb - 1 of data 
[311 ] (while this limit corresponds to a value of AR = 20, 
it does not depend strongly on AR) . Notice that we keep 
a uniform prior over AR for both choices of priors of 
the other two variables. We take a fiat prior over the 
SM nuisance parameters, whose value is however directly 
constrained by the likelihood - hence the choice of prior 
for those variables is unproblcmatic and it does not affect 
our results. 

The distribution of probability implied by our choice 
of priors for the MUED parameters and for some ob- 
servables is shown in Fig. [1] in one dimension, and in 
2-dimensional marginal distributions in Fig. [2j We ob- 
serve the expected uniform distribution in and AR 
for the flat prior choice in the left panels of Fig. [T] while 
the distribution on m/j is flat up to ~ 250 GeV and then 
it falls off sharply. This is a consequence of the fact that 
the radiative corrections to the KK Higgs are negative 
and proportional to the mass (m/j) of the SM Higgs 0]. 
The LEP limit on the Higgs mass is around 115 GeV, 
therefore leaving an allowed mass range 115-250 GeV. 
Therefore, for a given value of there is a value of 

the SM Higgs mass for which the charged KK Higgs be- 
comes lighter than the KK photon [32j , and it takes over 
the role of the KK photon as a DM candidate. As clearly 
charged DM is not allowed from cosmology, we discard 
points in which this happens. This sets an upper limit 
to the mass of the SM Higgs (to/, < 400 GeV, but with a 
very low probability above~ 250 GeV) as observed 
in Fig. [T] In the right-hand side panels of Fig. [TJ we can 
observe the impact of the log prior which disfavours large 
values of to/j and R . 

The distribution of the relic abundance under both pri- 
ors shows that very small values are not realized, i.e. the 
prior density goes to for SIdm^ 2 < 0.05. This is ex- 
plained by Fig. [21 where it is shown how the relic abun- 
dance is tightly correlated with which controls the 
level of degeneracy between the masses in the KK spec- 
trum. At lower values of R , the spectrum becomes 
more and more degenerate, hence annihilation is more ef- 
ficient and the relic abundance is reduced. However, since 
our prior includes a lower limit R^ 1 > 280 GeV due to 
the Tevatron constraints, as explained above, this leads 
to a lower limit in the distribution of ilj^^h 2 from the 
prior. While particle physics alone does not provide an 
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SM (nuisance) 


Mean value Standard deviation 


Ref. 


parameter 


/! G 




M t 


173.1 GeV 1.3 GeV 


[3JJ 


m b (m b ) MS 


4.20 GeV 0.07 GeV 


[38] 


a a {Mz) MS 


0.1176 0.002 


[38] 


l/a em (M z ) MS 


127.955 0.03 


[39] 



TABLE I: Experimental mean \i and standard deviation a 
adopted for the likelihood function for SM (nuisance) param- 
eters, assumed to be described by a Gaussian distribution. 





£l £2 £3 


£i 


5.78 x 10 6 -1.71 x 10 6 -4.65 x 10 6 


(■i 


1.39 x 10 6 8.93 x 10 5 


£3 


5.01 x 10 6 



TABLE II: EWPO covariance matrix employed in the analy- 




upper bound on R 1 , the thermal relic density of LKPs 
grows with and LKPs would overdose our universe 
for R~ x > 1.5 TcV Q. Finally, the prior distribution of 
both the spin-dependent and the spin-independent cross 
sections is relatively flat and spans several orders of mag- 
nitude. This range is to be compared with the much 
tighter range in the posterior (see below Fig. [3]), which 
means that the posterior distribution for those quanti- 
ties (to be discussed in detail below) is dominated by the 
likelihood. 

The likelihood function is constructed as follows. For 
each of the SM parameters, we assume a Gaussian likeli- 
hood with mean and standard deviation as given in Ta- 
ble |H To constrain the MUED parameters, we use data 
from electroweak precision observables (EWPO) which 
can be interpreted as constraints on the parameters given 
by the set [4fj| 

e = {ei,£2,£3}- (17) 

The maximum likelihood (ML) value of 9 obtained from 
LEP1 experiment data is [4l[ cml = {5 x 10~ 3 ,— 8 x 
10" 3 , 4.8 x 1CT 3 }. The likelihood function from EWPO is 
then modeled as a multi-dimensional Gaussian centered 
at the observed ML values, 

— 21ii£ewpo = Xewpo = ( e — e ML)*C~ (cml)) (18) 

where the covariance matrix C is given in Table HU 

We also include constraints from the WMAP 5-years 
measurement of the cosmological DM relic density [42j, 
which give for SIdm^ 2 a mean value /Avmap = 0.1099 and 
a standard deviation ctwmap = 0.0062 (notice that us- 
ing the updated WMAP 7-years values would not change 
our result considerably). When assuming that the LKP 
makes up the whole of the DM, we impose a Gaussian 
likelihood with the above mean and standard deviation, 
to which we add a 10% theoretical error in quadrature. 



We shall be interested in relaxing the requirement that 
all of the DM is made of LKPs, and in this case we use 
the WMAP measurement only as an upper bound. We 
show in the Appendix that in this case the correct effec- 
tive likelihood is given by the expression 

/>oo 

Aymap(^kk/i 2 ) = A) / e'^-^x^dx, 

(19) 

where Co is an irrelevant normalization constant, r* = 
/^wmap/cwmap an d Hkr^ 2 is the predicted relic den- 
sity of the LKP as a function of the MUED and SM 
parameters being considered. Notice that this is slightly 
different from what is usually adopted in the literature, 
namely either a sharp upper bound say 2cr above the 
WMAP mean, or a one-sided Gaussian which starts to 
drop at the WMAP mean and is fiat below (see Fig. QT]) . 

The total log-likelihood is thus given by the sum of the 
log-likelihoods defined above, i.e. 

— 2 In Aot = Xtot = Xewpo + Xsm + Xwmap- (20) 

The posterior distribution P(Q\d) is determined nu- 
merically by drawing samples from it. Markov Chain 
Monte Carlo techniques can be used to this aim, but 
in this paper we employ the MultiNest code, which im- 
plements the nested sampling algorithm (For a detailed 
description of the algorithm, see [H, [U, H^j). To per- 
form our statistical analysis, we use a modified version of 
the SuperBayeS code [3a, 13a] 1 which includes the Multi- 
Nest algorithm. Compared to standard MCMC methods, 
MultiNest provides a higher efficiency, guarantees a bet- 
ter exploration of degeneracies and multimodal posteri- 
ors and computes the Baycsian evidence as well (which is 
difficult to extract from MCMC methods). In our Multi- 
Nest scans, we use 20,000 live points and a tolerance fac- 
tor 0.5. We collect a total of about 220,000 samples from 
the posterior, which guarantees an adequate exploration 
of the parameter space. 



IV. RESULTS 

A. MUED parameter constraints 

We begin by showing in Fig. [3] the constraints on the 
MUED parameters and on some of the observable quan- 
tities. We plot the posterior distribution (blue), obtained 
by marginalizing over the posterior in the dimensions not 
show, and the profile likelihood (red), which maximizes 
the likelihood function over the variables not shown. The 
left panel is for the flat prior choice, while the right panel 
is for the log prior choice. It is clear that there is very 
little prior dependency and that the input parameters 



1 See www.superbayes.org 
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FIG. 1: Prior distributions for the input variables and some observables for flat priors on (mj, 7? 1 ) (left panel) and log priors 
(right panel), which are uniform in (log(m^), log(7? -1 )). 
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FIG. 2: Equal- weight samples from the prior for flat priors (top panels) and log priors (bottom panel). Density of points reflects 
prior probability density. 



are well constrained by the data. Therefore, from now 
on we will only show results from the flat prior choice. 
The profile likelihood is also in good agreement with the 
posterior, which signals that one expects little prior de- 
pendency. Hence our results can be deemed to be robust 
with respect to changes in the choice of priors and sta- 
tistical approach. 



The posterior distribution for peaks near the the 
best fit value R^ 1 = 641.6 GeV. The current Tevatron 
limit on as we have seen, is 280 GeV but by the 
end of 2011, Tevatron is expected to have 100 times more 
data, pushing up the limit closer to our best fit. By that 
time, the LHC should have collected 1 fb _1 of data, and 
it should therefore be able to discover MUED or at least 
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Parameter 


Mean 


Best fit 


68% range 


95% range 


LKP: the sole constituent of DM 


m h (GeV) 
R- 1 (GeV) 
AT? 


198.4 
640.9 
55 


215 
641.6 

38 


[173 : 222.3] 
[574.1 : 707.5] 
[23 : 86] 


[ 135.3 : 233.8] 
[536.5 : 843.5] 
[12 : 98] 


m 7 (GeV) 

^KK h 

\ogi4 1 ( P b)) 

log(a p SD (pb)) 


641 
0.115 
-11.1 
-5.7 


642 
0.111 
-11.2 
-5.7 


[574.7 : 707.8] 
[0.1 : 0.128] 
[-11.4 : -10.8] 
[-6 : -5.5] 


[537.3 : 843.4] 
[0.091 : 0.145] 
[-11.7 : -10.5] 
[-6.3 : -5.2] 


LKP: subdominant constituent of DM 


m h (GeV) 
R- 1 (GeV) 
AT? 


224 
602.9 

55 


226.7 
607.4 
66 


[202.4 : 245.4] 
[528.9 : 677.2] 
[25 : 86] 


[163.6 : 265.4] 
[477.1 : 795.5] 
[12 : 98] 


m 1 (GeV) 

i og (4^ ( P b)) 

log(a p SD C (pb)) 


603.5 
0.08 
-11.3 
-5.8 


607.9 
0.08 
-11.4 
-5.9 


[529.7 : 677.5] 
[0.057 : 0.108] 
[-11.6 : -11] 
[-6.1 : -5.5] 


[478 : 795.4] 
[0.035 : 0.127] 
[-11.8 : -10.6] 
[-6.3 : -5.2] 



TABLE III: Posterior mean and best fit values for the in- 
put MUED parameters and some relevant observables, both 
for the case where the LKP is the sole constituent of DM 
(top section) and where it is allowed to be a subdominant 
component (bottom section). We also give the 68% and 95% 
Bayesian equal-tails credibility intervals. While these figures 
are for the flat prior choice, the log prior choice gives very 
similar results and is therefore not shown. 



rule out the best fit (see the reach of the LHC in the 
right panels of Fig. 0]). As for the mass of the Higgs, the 
posterior peaks near the best fit value nih = 215 GeV, 
for which the Higgs dominantly decays into W + W~ and 
ZZ . This mass range of the Higgs is challenging for the 
7 TeV LHC with 1 fb" 1 . 

We do not find any constraints on the value of Ai? (see 
also Fig. 0|) . This can be understood as follows. In gen- 
eral, a change in Ai? modifies the mass spectrum, but the 
dependence is only logarithmic and it affects masses of 
strongly interacting particles only at the order of ~ 10% 
or less. For electroweak particles, the dependence is al- 
most flat in variation of AR. For instance, the KK lepton 
mass changes by 2% or so from AR = 10 to Ai? = 40 
and the KK photon is barely affected. In the computa- 
tion of the relic density, the dominant contribution arises 
from self annihilation of the KK photon and coannihila- 
tion with S , J7(2)v^-single KK leptons, while coannihila- 
tion with the KK quark is only subdominant (one rea- 
son is the coupling strength (hypcrchargc) and the other 
reason is the heaviness of the KK particles). Therefore 
the effect of a variation in Ai? is expected to be small 
in the relic density calculation. Similarly, the DD and 
ID are expected to be rather insensitive to variation in 
Ai?. Therefore, as far as the constraints considered here 
are concerned, only the values of i? _1 and mh are im- 
portant. However the collider phenomenology may be 
affected by the value of Ai?, and this will require a fur- 



ther, dedicated investigation. We notice that our results 
(and errorbars on the parameters) fully account for the 
lack of constraints on the value of Ai?. 

Fig. 0] shows 2D correlation plots for the MUED pa- 
rameters, both for the Bayesian posterior (top row, for 
the flat prior choice) and the profile likelihood (bottom 
row). Contour delimit regions of 68% and 95% prob- 
ability. We see also in this figure the relatively good 
constraints on i? _1 and rrih, and the lack of constraints 
on Ai?. We stress once more the reassuring agreement 
between the posterior and the profile likelihood, which 
implies little dependence of the priors. In Fig. [5] we plot 
the ensuing favoured regions for some of the observables. 

In Fig. [6] we show the mass spectrum of the first KK 
level for the best fit of the flat prior scan (the log prior 
case is very similar) , under the assumption that the LKP 
makes up the whole of the DM. The KK bosons (gauge 
(in green) and Higgs bosons (in magenta)) are shown in 
the left column, while the first two generations of quarks 
(in blue) and leptons (in red) are shown in the middle 
and the third generation in the right column. KK parti- 
cles denoted by lower (upper) case are singlets (doublets) 
under SU(2)w The mass spectrum and decay patterns 
from our best fit agree well with those shown in literature 
J3,B[S3] but here the corresponding scales are R" 1 = 642 
GeV and m h = 215 GeV. 

A natural question to ask is whether EWPO or relic 
density constraints influence more significantly the pos- 
terior. To understand this, we have performed two addi- 
tional analyses, keeping EWPO constraints and discard- 
ing the relic density bound in the first case, and vicev- 
ersa in the second. As expected, removing the bound on 
EWPO opens up the parameter space at small values of 
i? _1 , which results in a fairly flat pdf across the whole 
range of allowed Higgs mass values, and in pronounced 
volume effects on i? _1 , for which the posterior distribu- 
tion disagrees with the profile likelihood. Removing the 
constraint on the relic density completely (as opposed to 
setting an upper bound on this quantity, as discussed be- 
low), and keeping only EWPO constraints, allows larger 
values of R , which has the effect of pushing down the 
predictions for direct detection cross sections. 

To relax the strong assumption that the 71 particle 
makes up the whole of the relic DM, we adopt the upper 
limit on the relic abundance, represented by the likeli- 
hood in Eq. (fL9]) . This modifies the posterior in the f2 7 
vs m 7 plane in an obvious fashion, since we are allow- 
ing the case f2 7 < ft dm- The resulting constraints and 
corresponding favoured regions for the observables are 
shown in Figs. [7] and [5J In both those figures the EWPO 
constraints have been applied. There are small quanti- 
tative differences with respect to the case where the 71 
makes all of the DM, and the implications for the reach 
of the LHC and the mass of the Higgs are qualitatively 
unchanged. 

In terms of the best-fit x 2 , we find that the relic den- 
sity constraints contributes about 0.01 units to the total 
X 2 for the best fit point, which means that the WMAP 
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FIG. 3: Global constraints on the MUED parameters for two different choices of priors, assuming the LKP is the sole constituent 
of DM. The red cross gives the best fit, the vertical line the posterior mean. The horizontal blue/green bands give the 68%, 95% 
marginalized Bayesian posterior intervals; the red/green bands represent the 68%, 95% confidence intervals from the profile 
likelihood. There is only a very mild dependence in the constraints on the prior used or the choice of statistics. 



value can be reproduced very well by the model (both 
when it is taken as a Gaussian constraint and as an upper 
bound). The EWPO constraints contribute a x 2 ~ 1.6 
at the best-fit point. So our best-fit \ 2 is approximately 
1.6, for 1 nominal degree of freedom (4 data points for 3 
free parameters; we count 4 data points as the EWPO co- 
variance matrix has 3 independent eigenvalues, plus the 
relic density constraint). However, as we have mentioned 
above, the parameter AR is effectively unconstrained by 
the data, so it is not clear whether it should count in 
the computation of the number of degrees of freedom. 
In summary, our best-fit x 2 = 1.6 is statistically accept- 
able both when counting 1 degree of freedom in the fit, 
or (even more so) when discounting AR and therefore 
assuming 2 degrees of freedom. 



B. Prospects for MUED discovery 

We now move on to discuss the implications of our 
results for prospects of various experimental approaches 
to discover MUED. 

The rate of events in a direct detection experiment 
is obviously proportional to the product of the spin- 
independent cross-section, cr^ 1 , times the local density 
of the LKP, pj. This quantity can be smaller than /jdm, 
in the case where the LKP is not the only constituent 
of DM. In order to assess the prospects of detection, we 
therefore multiply cr^ 1 by £, i.e. the ratio between the 



local KK density and the local DM density, which, fol- 
lowing Ref. J44[, we assume to be equal to the ratio of 
the cosmic abundances of the two species, £ = p 7 //?dm = 
f2 7 /f2rjM- For 51dm we adopt the central value of the 
WMAP determination, while for pdm we adopt, follow- 
ing Ref. [ill, the value p x = 0.385 GeV cm -3 (see also 
[46H49j). We note that the actual DM local density is 
probably larger, due to the larger density of DM in the 
stellar disk [431, but we do not take this into account in 
order to be conservative. As one can see from Figs.[S]and 
[HI the 2-tr contours of the posterior lie below the sensi- 
tivity even of future experiments such as XenonlT [5fj| . 
An experiment attempting to probe KK DM in MUED 
should therefore be much bigger. The fact that the pos- 
terior is concentrated over one order of magnitude in tr^ 1 , 
and that it lies right below the sensitivity of 1 ton experi- 
ments, suggests that an eventual generation of 10-ton ex- 
periments, would be able to probe most of the favouted 
parameter space. 

The ID posterior of <Jp D , which is the key-quantity 
for indirect DM searches with neutrino telescopes, shown 
in Figs. [3] and [3 allows us to make a robust predic- 
tion on the prospects for detecting KK DM with the 
IceCube telescope, currently under construction at the 
South Pole, and already taking data. In fact, recent 
analyses of the sensitivity of IceCube to DM particles 
(see e.g. [5lj ) estimate the minimum cross-section that 
can be probed by this experiment to be k 5 x 10~ 5 pb or 
larger, in the relevant range of masses and after 5 years 
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FIG. 4: 2D global constraints on the MUED parameters for the case where the LKP is the sole constituent of DM. The top row 
shows 68% (yellow) and 95% (blue) regions from the posterior pdf (assuming fiat priors), while the bottom row gives confidence 
regions from the profile likelihood, with 68% confidence level region in yellow and 95% in red. We notice that the two statistics 
agree very well. The encircled cross gives the location of the best fit, the filled black dot of the posterior mean. In the mi, vs 
7? _1 plane we plot in magenta/green the 68% (inner contours) and 95% (outer contour) regions for the case where AR is fixed 
to 20 and 80, respectively. In the AR vs figure we show the LHC reach with 7 TeV and 1 fb _1 of integrated luminosity in 
the trilepton channel. 



of data taking. However, the posterior pdf for cr^ peaks 
one order of magnitude below this value, and it rapidly 
decreases for larger masses. Therefore we conclue that 
MUED searches at IceCube are unlikely to be successful. 

Turning now to the prospects at colliders, by the end of 
Run II, Tevatron is expected to deliver more than 10 fb _1 
of data and will greatly improve the current bound (280 
GeV at 95% C.L.), making it closer to our best fit point. 
From the LHC side, the reach for level 1 KK particles in 
MUED has been calculated in Q, where the gold-plated 
4£ $t signature is considered. The 4 leptons are obtained 
from the decay of KK Z, which is produced by the decay 
of KK quarks. This is quite similar to the production 
of the second lightest neutralino in supcrsymmetry. In 
MUED, however, the branching fraction of KK Z into 2 
leptons is large (1/6 for each generation) and the pro- 
duction cross sections of KK gluon and KK quarks are 
5-10 times larger than those in SUSY [H EM- The 14 
TeV LHC can probe MUED up to R' 1 ~ 1.5 TeV (1 
TeV) with 100 fb _1 (1 fb _1 ). A compactification scale 



of R^ 1 ~ 600 GeV (close to our best fit point) could be 
discovered or ruled out by the 14 TeV Run with 100 pb 
_1 H- The prospect for discovery of level 2 KK particles 
is discussed in [52|, [54| in terms of dilepton resonance. 
The reach is worse than the level 1 case due to the heav- 
iness of level 2 particles. Very recently, the reach at 7 
TeV LHC has been studied in J55|. It turns out that 
the opposite sign dilepton channel is the most promising 
discovery mode with 1 fb^ 1 of data. It is shown that 
MUED can be discovered if is less than 700 GeV, so 
this kind of search would be able to probe our la region 
for R- 1 . 



C. Distinguishing the MUED scenario from the 
CMSSM with direct detection 

We now turn to the question of how to distinguish a 
MUED scenario from a supersymmetric one, for which 
we will take the paradigmatic case of the constrained 
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FIG. 5: 2D correlations among some relevant observables and MUED parameters, for the scenario where the 71 particle is the 
sole constituent of the DM, with colour coding as in Fig. [4] We show the Bayesian posterior (top row) and the profile likelihood 
(bottom row). In the central and right-hand panels we display the reach of future direct detection experiments. 
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FIG. 6: Best fit mass spectrum of the first KK level from our 
global, assuming the LKP is the sole constituent of DM. 



minimal supersymmetric standard model (CMSSM) for 
simplicity. We briefly summarize in the following the 
approach taken here to constraining the parameters of 
the CMSSM (for details, see Ref. [561). 

Apart from the scalar mass mp, the gaugino mass 



and the trilincar coupling Aq assumed to be universal 
at Mqtjt, the CMSSM can be parameterized in terms 
of the bilinear scalar coupling B, the usual Higgs mass 
term in the superpotential and the SM-like parameters 
s. The latter include the 5(7(3) x 5(7(2) x (7(l)y gauge 
couplings, 53,3,3', and the Yukawa couplings, which in 
turn determine the fermion masses and mixing angles. 

In Ref. [56| was shown that considering M^ p as exper- 
imental data in the likelihood one can integrate out /1 via 
marginalization. This procedure automatically accounts 
for the fine-tunning in the sense that the posterior distri- 
bution is penalized in regions of the parameter space with 
large fine-tunning. Similarly the Yukawa couplings are 
easily integrated out when they are profitably traded by 
the physical fermion masses. Besides, it is highly advan- 
tageous to trade the initial B— parameter by the derived 
tan j3 parameter which is defined as the relative value of 
the two expectation values of the two Higsses. 

The resulting posterior in function of the usual vari- 
ables {mo, rrii/2-, Ag, tan /3)} introduces a global Jacobian 
factor in the posterior distribution which carries the pe- 
nalization of fine-tuned regions. Let us stress that the 
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FIG. 7: Global constraints on the MUED parameters for the flat prior, dropping the assumption that the 71 particle is the sole 
constituent of DM. The results with the log prior are very similar and are therefore not shown. 




FIG. 8: As in Fig. [4] but dropping the assumption that the 71 particle is the sole constituent of DM, and imposing only an 
upper bound instead. 
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FIG. 9: As in Fig. [5] but dropping the assumption that the 71 particle is the sole constituent of DM, and imposing only an 
upper bound instead. 



Jacobian is not "subjective" at all. Thus 

P(gi, m/,m , m^, A , tan p\ d) = 
Jl^nz P{gi,yf,m ,m 1 /2,A ,B,n = n z ) , (21) 

where J is the Jacobian of the transformation 
{fj,,yf,B} —> {Mz,mf, tan (3} (see Ref. for an 

explicit expression for J), \iz is the value of /i that re- 
produces the experimental value of Mz for the given val- 
ues of {s, mo, fnx/2, Aq, B} and P(s, m, M, A, B, /j,) is the 
prior in the initial parameters (For a detailed discussion 
on the chosen priors see Ref. 153). 

One of the main consequences of this approach is that 
the results exhibit a remarkable robustness under changes 
of the priors (see Ref. (57|), showing an absence of de- 
pendences on the initial chosen ranges for the CMSSM 
parameters. Moreover the results are compatible with 
likelihood based analyses [58j . 

Fig. [10] compares the favoured regions for the spin- 
dependent and spin-independent scattering cross section 
for the MUED and the CMSSM (see also Ref. |&|, where 
a similar analysis is performed). The experimental data 
used in constraining the latter are given in Table 2 of 
[5?}. Regions in light green (dark green) are within the 
reach of the LHC with 7 TeV and 1 fb" 1 (with 14 TeV 



and 100 fb 1 , respectively) for both models, while red re- 
gions are outside the reach of the LHC. Thus we can see 
that with 14 TeV and 100 fb" 1 the LHC is going to probe 
the whole of the favoured region for the MUED scenario. 
Also shown in Fig. [TU] are the sensitivities of various ex- 
isting and upcoming direct detection experiments. One 
sees from this plot that the detection of DM off spin- 
independent targets would point towards SUSY, rather 
than KK, DM, which is consistent with the findings of 
Ref. [l^ . The detection of KK DM in fact appears very 
problematic in astroparticle experiments. As we have 
seen in the previous section the spin-dependent coupling 
are such that the neutrino flux from DM annihilations in 
the Sun fall below the sensitivity of IceCube, even after 
5 years of data taking. 

The only viable search strategy appears to be the 
detection in an experiment sensitive to spin-dependent 
cross-section with a large exposure. The COUPP collab- 
oration, for instance, has been operating an ultraclean, 
room-temperature bubble chamber containing 1.5 kilo- 
grams of superheated CF3I, that produced interesting 
limits on the spin-dependent coupling [(32j . The plans for 
the future include the operation of a 60 kg chamber at 
Snolab, that could allow a substantial improvement in 
sensitivity. In Fig. [TU]we show for reference the reach of 
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FIG. 10: Favoured regions for the spin-independent and spin-dependent scattering cross section for the MUED scenario (bottom 
left cloud) and for the CMSSM. Points are equally-weighted posterior samples for each model, accounting for all relevant present- 
day constraints. Dark green (light gree n) regions are within the reach of the LHC with 7 TeV and 1 fb (with 14 TeV and 
100 fb -1 ) integrated luminosity [oOl J6lj[. while red points are outside the reach of the LHC. Closed black contours delimit the 
95% region for each model. Dashed lines give the approximate reach of future direct detection probes. 



the 60 kg version of COUPP, as well as the case of a 1 
ton chamber. If existing techniques for spin-dependent 
detection turn out to be scalable to such large volumes, 
and if one is not limited by some form of background, 
then KK DM could be within their reach. 



V. DISCUSSION AND CONCLUSIONS 

We have discussed the prospects for detecting KK 
DM at accelerators and with DM experiments with a 
a Bayesian analysis of the minimal UED scenario. Wc 
have derived in particular the most probable range of 
mass and scattering cross sections off nucleons, keep- 
ing into account cosmological and electroweak precision 
constraints. The value of the three free parameters of 
the model at our best fit point are Rr 1 = 641.6 GeV, 
AR = 38 and m h = 215 GeV if the KK DM explains 
all of the DM in the universe and i?" 1 = 607.4 GeV, 
AR = 66 and m h = 226.7 GeV if KK DM is assumed 
to be a subdominant constituent. As we have seen, the 
current Tevatron limit on R~ 1 is ~ 28QGeV, but by the 
end of 2011, Tevatron is expected to have 100 times more 
data, pushing up the limit closer to our best fit. By the 
time, the LHC should have collected 1 fb -1 of data, and 
it should therefore be able to discover MUED or at least 
rule out the best fit. For the two DM scenarios our best 
fit points for the Higgs mass are ~ 215 — 227 GeV, for 
which the Higgs dominantly decays into W + W~ and ZZ. 
This mass range of the Higgs is challenging for the 7 TeV 



LHC with 1 ftr 1 . 

Our analysis has dramatic consequences for the de- 
tcctability of the MUED scenario with astrophysical DM 
experiments. Figs. [5] and |H] clearly show that the 2-sigma 
contours in the a^ 1 vs. mass plane fall below 10~ 10 pb, 
i.e. even beyond the reach of future ton-scale experi- 
ments. This implies that if new particles are actually 
found with direct detection experiments, they are un- 
likely to be associated with KK DM. Direct detection is 
however not hopeless, provided that current experiments 
with to spin-dependent targets rapidly improve their sen- 
sitivity. We have seen that experiments such as COUPP 
might probe the relevant portion of the parameter space 
if they can go beyond the upcoming scale of 60 kg. The 
analysis presented here therefore provides an additional 
motivation to build such detectors, in case the MUED 
scenario is discovered at accelerators, in which case one 
could perform a combined analysis of accelerator and di- 
rect detection data, following the approach suggested in 
Ref.fUl. 



Indirect detection prospects are not very promising, 
with the most probable flux from of neutrinos from KK 
DM annihilations at the center of the Sun below the 
sensitivity of IceCube, even after 5 years of observa- 
tion. We haven't discussed explicitly the possibility of 
detecting gamma-rays, anti- matter or synchrotron emis- 
sion from KK DM annihilations in the halo, because the 
prospects for detection depend strongly on the assump- 
tions made on astrophysical parameters, and when con- 
servative choices are made for these parameters, the pre- 
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dieted fluxes are below the astrophysical backgrounds. 
This is easy to understand, since the annihilation fluxes 
typically scale like ~ av jmt , and the annihilation cross 



section is av ~ m 2 , it follows that 



, 4 . Therefore, 



» 7 , _~ — r ,.„ 7 

given that the most probable range of mass is centered 
around the relatively large value of 600 GeV, all annihi- 
lation fluxes are quite suppressed [III [64| . 

Although we discussed a minimal version of the UED 
scenario, various extensions beyond MUED have been 
suggested and their rich phenomenology of Kaluza-Klein 
dark matter has been investigated in Rcfs. 
[651 - 4671 ] . An analysis of these non- minimal scenarios will 
be the subject of a dedicated forthcoming paper. 
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Appendix A: Derivation of the relic density upper 
bound likelihood 



In this appendix, we derive the likelihood given in Eq. 
(TT9"|) . There exist in the literature various expressions 
for the likelihood function in the case where the WMAP 
result is taken to be only an upper bound to the DM 
density. Often, the likelihood is taken to be flat up to 
an arbitrary cutoff value (e.g., the 95% upper r ang e of 
the WMAP likelihood) and zero above it. Ref. [631 a d" 
vocated using a likelihood function which is flat below 
the WMAP central value, and falls off as a half- Gaussian 
above it. Here we derive the correct expression, which is 
functionally slightly different from what has been previ- 
ously used. 

We define the following shortcut notation: ojkk = 
Okk/i 2 is the relic density of KK particles, while wdm = 
floM^ 2 is the relic density of all dark matter, which might 
comprise a secondary component beside the LKP, i.e. 
wkk < wdm- The WMAP measured mean value is given 
by /xwmap, and its uncertainty is ctwmap- We thus want 
to determine the effective likelihood 



^WMAp(wkk) = p(a*WMAp|wkk) — / p(A i WMAp|wDM)p(^DM|'^KK)dwDM 



(Al) 



r 



where p^wmapIwdm) is a Gaussian in u>bm with mean 
A*wmap and standard deviation <t\vmap, i-e. wdm ~ 
A/"(/i W MAP, CTwmap)- In order to determine p(wdm|wkk), 
we use Bayes Theorem to obtain 



p(w D m|w K k) =J3(w K k|w D m) 



p("dm) 
p(wkk) 



(A2) 



On the RHS of Eq. (|A2j) the first term is the conditional 
probability for the LKP relic density given a specified 
total DM density. Since we are considering the case 
wkk < udm and nothing else is known about the rel- 
ative densities between the LKP and a secondary dark 
matter component, we set 



p(wkk|^dm) 



otherwise. 



(A3) 



To specify the priors p(wdm),p(wkk) in Eq. (|A2|) . we ap- 
peal to the principle of indifference. Lacking any other 
information about the relative densities of the LKP par- 
ticle and the total DM density, we should take the two 
priors to be equal for both components. The ratio there- 
fore must be 1 everywhere except in the unphysical region 
of negative energy density, where we set it to to enforce 
positivity of the energy densities. Thus we have 



p(^dm) _ J 1 if < wkk.wdm < 0, 
p(wkk) I otherwise, 



(A4) 



where f2 is some large cut-off value whose precise value 
is irrelevant for the end result, as it will be shown below. 
Using Eqs. (|A2HA4p into Eq. JSTJ we obtain, taking the 
limit f2 — > oo: 



14 



£wmap(wkk) = lim —— 

o^oo v27ra"wMAP 



exp 



1 (^DM — MWMAp) 2 
2 



exp 



WMAP ^ w kk/<twmap 



1 



WMAP 
2 



W D J,0(tdKK — WDM)dwDM, 



(x - r*) 



r 



"Ma;, 



(A5) 
(A6) 
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which is Eq. (fT9|) (notice that although this expression 
is not normalized this is immaterial as we only need the 
likelihood up to an overall normalization constant). This 
effective likelihood is plotted in Fig. [11] where it is com- 
pared to the Gaussian likelihood for the case when the 
DM is made entirely of LKP. The effective likelihood is 
flat for lukk -C ^wmap, then falls off exponentially for 
wkk 3> Mwmap, as one would expect. Notice that the 
likelihood at wkk = MWMAP is precisely half is asymp- 
totic value for lokk "C /Avmap, which reflects the fact 
that for the WMAP central value we are agnostic as to 
which fraction of DM is made of KK particles. 



FIG. 11: Comparison between the WMAP likelihood when 
the LKP is the sole constituent of dark matter (green, Gaus- 
sian shape) and when the LKP is a subdominant component 
(red, upper bound). Both likelihoods are normalized to their 
peak value. 
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